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ABSTRACT 

We present deep two-colour photometry of two rich clusters at z — 0.18, A665 and 
A1689. We use these data to construct number counts as a function of magnitude 
in the two fields. By combining these counts with similar observations from a large 
area field survey we subtract the field contamination statistically to produce luminosity 
functions for the two clusters. Great care has been taken to achieve agreement between 
the photometry of these two samples. The cluster data are complete to a limiting 
magnitude of / = 22.5 or an absolute magnitude in the cluster of / = —18.0 {M* + 5). 
The luminosity functions of both clusters are well described by a Gaussian function for 
the bright galaxies, combined with a Schechter function at the faint end, similar to that 
required to fit the luminosity function in local clusters. The slope at the faint end of 
the Schechter function in both clusters is extremely steep in y, a — 2. A shallower 
slope is seen to the limit of the / data, indicating that the cluster population is rapidly 
blueing as we reach fainter. The excellent agreement between the form of the luminosity 
function i n our two distant cl usters, as well as agreement with the luminosity function 



given by Driver et al. (1994) for a single z = 0.21 cluster, indicates that this faint 
blue population is a general constituent of distant clusters. We compare our results 
with those from studies of local clusters. Depending upon the degree of fading (or 
disruption) of these faint blue galaxies, we tentatively identify their remnants with 
the low surface brightness dwarf galaxies which are the dominant population in local 
clusters. We discuss the possible role of this population as the source of most of the 
X-ray gas in rich clusters. 

Key words: cosmology: observations - clusters: individual: A665, A1689 - galaxy 
evolution - galaxies: photometry - galaxies: luminosity function. 



1 INTRODUCTION 

The form of the luminosity function (LF) of ga laxies is a 



critic al constraint on models of galaxy formation (Cole et al 



1994). W6 taeaiOita aCCUi'atCly thti fol'm of thU galaxy luttii- 



nosity mncLiun Lu fainl maguiLudcs icquiies Lhc delcimina- 
tion of distances to large numbers of faint galaxies. Observa- 
tionally this is a very expensive procedure to undertake for 
field galaxies, as the individual distances must be measured 
spectroscopically. However, by studying a rich cluster where 
the bulk of the population is at a single distance, the practi- 
calities of LF estimation are optimised. Unfortunately, this 
advanti.ge comes at a price; rich clusters are regions where 



we might expect environmental effects on the form of the 
luminosity function to have their greatest impact, evidence 
for which may alrea dy exist in th e observed morphology- 
density rela tionship (Dressier 1980; Whitmore, Gilmore, & 
Jones 1993). Nevertheless, by studying a number of clusters 
with a range of properties it may be possible to gauge the 
effect of the environment on their galaxy populations (itself 
an interesting issue) and hence determine the likely similar- 
ity of field and cluster LFs. It is for this reason that much 
recent work has been concentrated on rich clusters. 



A photographic study of 14 rich clusters by CoUess 

(1987) and another 9 clusters by Lugger (1986) concluded 
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that all clusters could be fit by a universal LF with a ^ 
-1.25 to Mv ~ -I9Q More recent work has come to some- 
what different conclusions. Kashikawa et al. (1995) carried 



out a photometric survey of both morphologically-classified 
subsamples and the composite LF's of four nearby clusters, 
including Coma. They found a faint end upturn which led 
to a rejection of a single Schechter function with a ~ —1.25 
as an adequate representation of the data; a combination of 
a Gaussian at the bright end and a Schechter function at 
the faint end being p referred. A faint end upturn in Coma 
was origi nally scon by lAbcU (19771) an d later [Metcalfe (19831 ) 
(See also Thompson fc Gregory (1993 ), Biviano et al. (199 



and Seeker & Harris (1996), all of whom also prefer the com- 
bined Gaussian-fSchec hter function descripti on of the LF). 
However, we note that Bernstein et al. (1996) found no up- 
turn in their counts in th e core region of Coma. 



In 



(1985^ 



a seminal paper Biiiggeli, Sandage, & Tammami 



mblihhed the LF derived fi 



-tire 



extensive studji 



of the galaxy population in the Virgo region. The advantage 
of working in the Virgo cluster is that its relative proximity 
means that it is possible to reach far down the LF. Using 
deep plate material they classified the galaxy populations 
morphologically and then derived luminosity functions for 
various types. The brighter cluster ellipticals apparently fol- 
low a Gaussian distribution, while the faint end of the LF 
is dominated by a population of dwarf ellipticals which im- 
part a steeply rising faint end slope t o the LF (a ~ —1.35) 



Later, [mpey, Bothun, fc Malin (1986 ) (see also Bothun, Im- 



pey, fc Vlalin 1991) surveyed both Virgo and Fornax for low 
surface bri ghtness (LSB) object s and found many galaxies 



missed by Binggeli et al 
to 



(198! 



steepening a yet further 
1.7 in Virgo (brighter than Mv ~ -12). While the 



breakdown of the LF into its separate morphological con- 
stituents was extremely illuminating, this type of study can 
only be pe rformed for nearby clusters from the ground (al- 
though see small et al. (1996) for a study of distant clusters 
using HST) and suffers from uncertain corrections for field 
contamination. 

A general picture is thus emerging of the form of the 
luminosity function for local clusters. The LF appears to 
be best described by combining a Gaussian distribution for 
the bright (giant) galaxies with a Schechter function for the 
fainter (dwarf) galaxies. The very steep slope of the faint 
component results in it being the dominant galactic com- 
ponent in the clusters. The evolution of this population is 
thus critical for understanding, for example, how the blue 
star forming galaxies seen in the cores of distant (z > 0.2) 



But cher -Oemler clusters (Butcher fc Oemler 1984 



Oemler 



1992) hive evolved by the present-day. These galaxies should 
still populate such regions, although their appearance will 
depend sensitively upon their subsequent star formation his- 
tories. One recently proposed mechanism for tran sforming 
them is "galaxy harassment" ( Moore et al. 199^ ). In this 
model star forming field galaxies falling into the cluster po- 
tential have their dark halos stripped, as well as undergoing 
tidal encounters with cluster galaxies. These processes re- 



* We adopt qo = 0.5 and Ho =50 km sec ^ Mpc ^. With these 
parameters 1 arcsec is equivalent to 3.91 kpc at 2 = 0.18. 



move the outer baryonic material from the galaxy disk and 
concentrate its remaining gas in a nuclear starburst. The end 
product of the process is suggested to be a dwarf spheroidal 
galaxy. Within the framework of such a model we would ex- 
pect the dwarf/giant spheroid ratio of the cluster to increase 
rapidly with cosmic epoch, as increasing numbers of infalling 
galaxies are transformed into dwarfs. 

Another area of research which is sensitive to the form 
of the faint end of the LF in distant clusters is the source 



of intracluster gas and its enrichment. Trentham (1994) has 
suggested that a large fraction of the hot intracluster gas 
seen in Coma may have originated in dwarf galaxies which, 
after an initial burst of star formation, expelled their gas into 
the cluster potential and subsequently faded to become the 
LSB dwarf galaxy population observed today (or disr upted 
completely) . To explain all the gas observed in Coma 



Tren- 



tham concluded that an a ~ —1.8 faint end slope to the LF 



is required in the distant clusters, moreover this steep popu- 
lation should be predominantly star forming dwarf galaxies. 
It is therefore of considerable interest to determine observa- 
tionally the form of the faint end of the galaxy luminosity 
function in distant clusters. 



Driver et al. (1994)[DPDMD] made the first attempt to 
measure the faint end of the LF in a distant cluster {z > 0.1). 
Their study showed evidence for a steep slope t o the LF in 
the range -I- 2 to M* -f 6. [Driver et al. (199^ ) obtained a 
total integration time of 2.4 ksec in B and i? of a 6 x 4 ar- 
cmin region in the distant cluster A963 {z = 0.206) using the 
Hitchhiker parallel camera on the 4.2-m WHT. They used 
the data to investigate the number counts in this region and, 
by correcting statistically for the field contamination from 
counts obtained with the same instrument in a number of 
blank fields, they attempted to measure the cluster galaxy 
LF down to i? ~ 24 (equivalent to M* + 6) . They found a 
steeply rising faint end slope, with a ~ —1.8 in a double 
Schechter parametrisation. The slope of the faint end of the 
LF, {dlog N /dm) ~ 0.3, was similar to the observed slope 
of the field counts. Thus their result could be reproduced 
by a simple zero-point shift between the magnitude scales 
of their field and cluster images, albeit of rather large am- 
plitude, Sm ~ 0.2. Unfortunately, owing to its operational 
mode the Hitchhiker system was only calibrated once a year, 
instrumental consistency has to be assumed, and thus it is 
difficult to gauge the likelihood of such an error (They quote 
an rms scatter on the zero point of ±0.1 mag). To address 
this issue DPDMD compared their photometry with data 
from a shallow photographic survey which appeared to rule 
out any large offset. 

To provide an independent check of the result reported 
by DPDMD we have analysed deep two colour photometry 
of two distant z — 0.18 clusters, using a similar technique 
to measure the cluster galaxy luminosity function. Our clus- 
ter dataset is similar in quality and depth to DPDMD 's, 
although roughly a factor of six larger in area for each clus- 
ter than that of DPDMD. To provide field counts with the 
same selection criteria and conditions as our cluster images, 
rather than relying on published counts from the literature, 
we have analysed a wide-field V and I imaging survey of 
blank fields. In addition we have taken great care to ensure 
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Figure 1. An / image of the rich cluster A665. The field of view 
is 12.3 X 12.3 arcmin (2.9 X 2.9 Mpc). North is top and East is 
left. Overlaid is a contour map of the X-ray emission from the 
cluster. The X-ray emission peaks on the central cluster galaxy 
and shows an asymmetric distribution with a plume off to the 
north-west. This structure is also visible in the distribution of 
red cluster galaxies. 
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Figure 2. The I image of the cluster A1689. The field of view is 
11.3 X 10.6 arcmin (2.6 X 2.5 Mpc) and North is top with East 
to the left. Again we have overlaid the cluster X-ray distribution 
as a contour map. The highly symmetrical X-ray peak coincides 
with the dense clump of bright galaxies which defines the cluster 
centre. 



the homogeneity of the photometric systems used for both 
the cluster and comparison blank fields. 

We discuss our observational dataset in Section ^ 
present our analysis in Section ^, discuss our results in Sec- 
tion ^ and give our main conclusions in Section |^. 



images taken from the ROSAT archive to illustrate the gen- 
eral morphologies of the clusters. The two clusters are both 
strong X-ray sources with luminosities in the 2-10 keV band 
of Lx = 1.2 X lO'"^ erg s sec~^ for A665 and L x = 2.0 x lO"^ 
ergs sec~^ for A1689 (Soltan & Henry 1983). In addition, 
both clusters have high velocity dis persions: ad — 120 1^1^^? 
km sec~^ (33 members) for A665 ( [Oegerle et al. 1991) and 
aci = 1848 ± 166 km sec~^ (68 meinbers) for A1689 (IGudc-' 



hus 198£). Finally, both clusters sh ow evidence of strongly 



2 OBSERVATIONS AND REDUCTIONS 



Before 



lensed features in their core regions ( Tyson, Valdes, fc Wenk 



discussing uui ubsyivaLiuus of Lhy clusLyis and Lhy 



199C; Kaiser et al. 1994; Tyson fc Fischer 1995) 



field we first describe the general characteristics of the two 
cluster^ observed. 



The characteristics of these clusters are very similar 
to the cD-dominated cluster A963 studied by Driver et al 



(1994). A963 is a richness class 3 cluster at z = 0.206 which 



2.1 



Tusfer'Tropernes" 



: 0.95 x lO*'^ ergs sec ~\ feoltan 
fc Henry 1983) and also contains two giant arcs (Lavery & 



is a strong X-ra y source {L^ 



The two clusters analysed in this study were A665 and 
A1689. These clusters both lie at 2 ~ 0.18 (A665 at z = 
0.182 and A1689 a.t z = 0.181). They are both optically 
rich with Abell richnesses of 5 and 4 respectively. Using 
our imaging data we estimate the blue fractions {fs from 
Butcher fc Oemler 1984) of both the clusters to be low, 
/s = -0.11 ± 0.19 for A665 and /s = 0.12 ± 0.05 for 
A1689 (c.f. /s = 0.09 ± 0.03 for A1689 from Butcher & 
Oemler 1984). They have somewhat different morphologies, 
with A665 having a single cD galaxy and A1689 containing 
a very dense central group of galaxies (See Figs |l| and |^. 
In the figures we have overlaid the ROSAT PSPC X-ray 



Henry 1985 ; Ellis et al. 1991). The cluster's blue fraction has 
been estimated as /s ~ 0.19 ± 0.05 by Butcher & Oemler 
(1984). 



2.2 Data Acquisition 

The V and I imaging of the two clusters used in our analy- 
sis was collected with the Prime Focus imager on the 2.5-m 
Isaac Newton Telescope (INT), La Palma. The comparison 
V and I observations of blank fields needed to correct for 
field contamination come from deep imaging with the f/1 
camera on the 3.9- m Anglo- Australian Teles cope (AAT), 

for more de- 



Siding Springs (See Lidman & Peterson 1996 
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tails) . In addition we obtained additional photometric imag- 
ing to tie the photometric systems of these two data sets 
together securely. These data were collected with the Prime 
Focus imager on the 4.2-m William Herschel Telsescope 
(WHT), La Palma. 

The imaging data on the two clusters, A665 and A1689, 
was collected using a 2k^ thick Ford CCD on the nights of 
22-28 February 1993. Due to a technical problem with the 
data acquisition system running the Ford chip a thick Ik^ 
EEV CCD was used on the night of the 23-24 February 1993. 
The seeing during this run varied between 1.4-2.3 arcsec, far 
in excess of the limit required for the intended programme 
of lensing observations. Accordingly, the current project was 
executed as a backup programme. This involved imaging in 
V and I of two clusters ai z ^ 0.18 to study their LFs over 
a wide-field. 

Our observing technique was to take multiple exposures 
(each of ~ 1000s) of the clusters, dithering the telescope 
by ~ 15 arcseconds between exposures. This reduced flat- 
fielding errors and allowed us to create a master fiatfield 
as explain ed in the next section. Observations of standard 
stars from Landolt (1992) were interspersed with the science 



Burstein fc Heiles (1984) 



observations throughout each night. The standards were ob- 
served across a large range in airmass to provide atmospheric 
extinctions and zero points for the various clusters and pass- 
bands. Furthermore, the colour range chosen for the stan- 
dard stars spans the same range as expected for the cluster 
galaxies allowing us to determine colour-terms for the vari- 
ous detector and filter combinations. Two nights during the 
run were photometric and calibration of all passbands was 
made using these nights. We also corrected our magnitudes 
for galactic absorption (reddenin g) using Ai = 0. 50 Ab and 
Av ~ 0.75Ab with Ab given by 
(see Table |). 



2.3 Data Reduction 

The reduction of our data followed standard procedures for 
the analysis of deep imaging. In particular, we chose to cre- 
ate super-flats using the data frames themselves to allow us 
to best match the characteristics of the flatfields to the sci- 
ence images. The reduction was all performed with standard 
IRAF routines and the procedure in detail was: 

(i) The frames were bias subtracted and trimmed. The 
median bias level was obtained from the overscan region of 
the chip and subtracted off. The images were then trimmed 
to remove the overscan strips. 

(ii) Initial fla tfielding was perfo rmed using twilight fiats, 
(ui) FOCAS dValdes et al. 1983|) was used to detect bright 



objects in the frames. These were then removed and replaced 
by sky values drawn from regions around the objects. 

(iv) For a given frame, all the other cleaned frames in that 
passband taken on the same night were median combined to 
create a super-fiat. 

(v) The initially fiatfielded frames were then fiatfielded 
using the super-fiat. 

(vi) The images were aligned and geometrically re- 
mapped using sub-pixel sampling to a single basis frame 
(centred on the cluster), and then combined using a clipped 



average algorithm in the IRAF task IMCOMBINE to give 
final images for each passband/night. 

(vii) All final exposures from the different nights, in a 
given passband, were co-added to provide a single frame of 
each cluster in each passband. 

The characteristics of the final datasets are summarised 
in Table ^. Our cluster data covers a region of approximately 
12 X 12 arcmin or 2.8 x 2.8 Mpc in each cluster to a depth 
of / ~ 22.5 and V ~ 24.0. These are both equivalent to 
M* -I- 5 at the cluster redshift. The average seeing ranged 
from 1.7 to 2.1 arcsec FWHM (Table |l|). We use the mean 
colour of the elliptical galaxy population ((V^ — /) ~ 1.5) 
to determine the colour correction to transform our magni- 
tudes onto the Landolt system. We include a contribution 
to the final magnitude errors from differences between the 
adopted colour correction and the true value arising from 
the observed range of galaxy colours in our fields. The final 
zero point errors were AI = 0.05 and AV = 0.05, obtained 
from a combination of extinction, colour and frame-to-frame 



2.4 Cluster Galaxy Catalogues 

Having acquired and reduced our images we next needed to 
analyse them to provide catalogues of object positions, mag- 
nitudes and colours. To achieve this we used the SExtractor 



analysis package (Berlin & Arnouts 199(3). This package de- 
tects objects using an isophotal threshold above the local sky 
(after convolution with a filter function), cataloguing those 
with areas exceeding a minimum value. SExtractor then de- 
blends the objects and produces a catalogue of their prop- 
erties. The object catalogue created includes information on 
the object positions, shapes, profiles and magnitudes. Tests 
to find the best detection parameters for SExtractor were 
run on small sections of both the V and I images of the 
final cluster frames. The final values we adopted are shown 
in Table ^ While we used SExtractor to detect objects and 
determine their centroids, to compare the number counts in 
our various fields we have chosen to use large-diameter aper- 
ture photometry (6 arcsec diameter or 23.5 kpc at the clus- 
ter distance). This photometry was performed with IRAF's 
PHOT package and we discuss it in more detail later. 

We corrected the catalogues for incompleteness a t the 
faint end by means of simulations (c.f. Bmail et al. 1995). We 
artificially generated a "typical" galaxy for each of our four 
cluster fields by co-adding a large number of bright galaxies 
(7 ~ 20). We then scaled the fiux in this "typical" object to 
obtain simulated galaxies of different magnitudes. We added 
200 such galaxies of a fixed magnitude at a time to each of 
our science frames. We then re-ran SExtractor on the re- 
sulting frames and noted what percentage of the additional 
simulated galaxies were being missed. The simulations thus 
include the effects of incompleteness due to both the faint- 
ness of the galaxies and also merging with brighter objects. 
Note that implicit in these simulations is the assumption 
that, over our range of interest of 1 — 3 magnitudes, faint 
galaxies have similar scalelengths to their brighter counter- 
parts. Even if this were not the case, and fainter galaxies 
were intrinsically smaller, the high values of seeing (~ 2 
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arcsec) encountered during these observations make it likely 
that all objects of differing magnitudes would be smoothed 
to a similar size. Table ^ gives the 50 per cent completeness 
limits of the catalogues obtained from these simulations. 

To correct the catalogues for stellar contamination we 
have used the relation between image si ze and brightness 



to separate the stars and galaxies (c.f. Lidman & Peter- 
son 19!|^. In a plot of isophotal radius versus magnitude 
two distinct populations are evident at bright magnitudes, 
with the more compact objects at any magnitude being the 
stars. For each frame a line separating the populations was 
determined visually and the objects more compact than this 
limit were classified as stars and removed from the counts. 
At faint magnitudes (/ > 18) the two populations begin to 
merge and it is more difficult to identify the stars. Hence it 
seems likely that some contamination will remain, although 
the fraction of stellar interlopers is small. We are primar- 
ily interested in the faint end of the LF and owing to the 
relatively shallow star number counts compared to galaxies 
the star contamination becomes significantly reduced as th e 
magnitude increases (see fig. 4 of Lidman & Peterson 1996). 



After star removal, the A665 and A1689 / catalogues and 
A665 and A1689 V catalogues contained 1689, 2035, 1659 
and 1615 objects above their respective 50 per cent com- 
pleteness limits. 

Before discussing the field counts used to correct our 
cluster observations we first discuss whether the SExtractor 
detection parameters that we used might cause us to miss a 
population of low surface brightness galaxies. 



2.5 Are Low Surface Brightness Galaxies Being 
Missed? 

CCD frames are less susceptible to missing low s urface 



brightness objec ts t han photographic p lates (see e.g. Bern 



3tein et al. 1996 



and Turner et al. 1993 



for more discussion) . 
Despite this, and to test the sensitivity of our object detec- 
tion algorithm to the adopted surface brightness limit, we 
also catalogued our cluster frames using a procedure simi- 
lar to that used by DPDMD in their analysis of A963. We 
selected a very low surface brightness limit for our object 
detection with the understanding that the ensuing sample 
would be strongly contaminated by noise objects. Rather 
than making an arbitrary choice for the detection threshold 
we chose to vary its value in the range 1.0-2.0 a and se- 
lected the level which produced the highest number of new 
detections. This is not simply the lowest threshold, owing 
to the effects of the deblending algorithm on merged galaxy 
images. 

To determine the excess number of real objects not de- 
tected by our standard procedure we first removed all those 
objects from the low-threshold (LSB) / catalogues which oc- 
curred in the standard versions of the catalogues. Next we 
used the fact that the reddest galaxies expected at the clus- 
ter redshifts {z = 0.18) will correspond to the colours of the 
spheroidal populations (E/SO) of the cluster {{V — I) ~ 1.6). 
Thus any real cluster objects will have colours bluer than 
{V — I) < 1.6 measured from the seeing-matched V and 
I images, or for our limiting magnitudes of / = 22.0 and 
/ = 22.5 they would have V < 23.5 and 24.0 respectively. 



The number of objects from the LSB catalogues which 
are brighter than I = 22.0 (22.5), are undetected using the 
standard algorithm, and yet have V magnitudes brighter 
than V = 23.5 (24.0) was 50 in A665 (2.5 per cent of the 
whole population) and 33 in A1689 (1.4 per cent) respec- 
tively. These proportions make a negligible difference to the 
cluster LF discussed later and so we conclude that our stan- 
dard object algorithm is not significantly biasing the object 
catalogue against low surface brightness cluster members, in 
so far as they are detectable in our data. 

2.6 AAT f/1 Field Imaging 

The blank field observations necessary to field-correct our 
cluster counts come from a wide-field V and / CCD sur- 
vey of the equatorial LOSS fields at lO'^"" and 13*"'=. The 
detector was a Ik^ Thomson CCD with 0.98 arcsec/pixel 
sampling, providing a large field of view (~ 17 x 17 arcmin). 
The data were all taken in photometric conditions and in 
seeing comparable to that encountered for our cluster ob- 
servations. These data were kindly supplied and reduced by 
Dr. C. Lidman. Dr. Lidman also provided internal calibra- 
tions between the various fields with estimated internal zero 



point errors of A ~ 0.03 (Lidman & Peterson 1996) 



2.7 



Calibration of the Comparison Field 
Photometry 



This final imaging dataset was acquired to tie the AAT and 
INT frames rigorously to a single photometric system. While 
both catalogs claim to be on the Landolt (1992), it appeared 
prudent to us to actually test this claim. For this purpose 
the INT cluster fields and a selection of the AAT blank 
fields were re-observed with a Ik^ thinned Tek-2 chip on the 
Prime Focus Imaging Facility of the WHT. The conditions 
were photometric but the seeing wa s very poor (~ 4 arc- 
sec). Observations of Landolt (1992) standard stars across 



a wide range in airmasses were also taken to provide robust 
photometric transformations. The science frames were re- 
duced in a standard manner using twilight fiatfields before 
being rebinned to either the INT (for cluster) or AAT (for 
field) pixel scales. The INT and AAT data were degraded to 
the same seeing as the WHT and 10 arcsec diameter aper- 
ture magnitudes then measured from galaxies in the fields. 
We determine magnitude offsets for the AAT data relative 
to the INT cluster fields. The AAT field data was deemed 
to have magnitude offsets of A{Iaat — I) = 0.07 ± 0.05 
and A{Vaat - V) = Q.U ± 0.07 relative to the INT sys- 
tem. The relatively large offset uncovered in the V band is 
disappointing given the claimed accuracy of the two cata- 
logs, but we reiterate that it was to check for exactly this 
sort of problem that we undertook these observations. The 
source of this offset has proved difficult to track down. One 
check on the relative magnitude scales in the INT cluster 
data is available from the colours of the E/SO sequence, this 
has a mean colour of {V ~ I) — 1.51 at V = 17.2 with a 
scatter of S{V — I) = 0.05 in both clusters. This is in rea- 
sonable agreement with the no evolution predicted colour of 
{V — I) = 1.57, indicating that the INT photometric is likely 
to have offsets of < 0.06. 



6 Wilson et al. 



We applied corrections to the field catalogues to put 
the field and cluster photometry all onto the INT system. 
By performing this calibration of our control fields we were 
able to ensure that the relative zero point errors between our 
cluster and field photometry were reduced to a minimum 
(~ 5 per cent random errors). 



2.8 Field Catalogues 



We analysed the AAT imaging data as in Section 2.4, with 



equivalent selection criteria. In our analysis we chose to use 9 
I fields (with seeing less than 2.1 arcsec) and 4 V fields (with 
seeing less than 2.6 arcsec). We employed similar simulations 
to those described in Section 2.4 to correct for incomplete- 



ness, this time generating a typical galaxy for each passband 
(the observing conditions and, in particular, the seeing for 
all the data frames taken with each filter was approximately 
similar so it was only deemed necessary to perform incom- 
pleteness simulations using one representative frame in each 
passband). The AAT field survey provided adequate statis- 
tics for number counts as faint as 7 ~ 22 and V ~ 22, close 
to the limiting magnitudes of our cluster datasets. Again, 
we isolated and then removed the stars by plotting isopho- 
tal radius versus magnitude. 

To measure the magnitudes of the objects on the clus- 
ter and comparison fields we chose to use simple aperture 
magnitudes. To provide reasonable estimates of "total" mag- 
nitudes these were measured in 6 arcsec diameter apertures 
centred on the object positions provided in the SExtractor 
catalogues. The background sky was measured in a wide 
surrounding annulus. To remove differential aperture cor- 
rections between the datasets we convolved all the images 
in a particular passband to the seeing of the worst, using a 
Gaussian filter, before measuring aperture photometry with 
IRAF's PHOT package. 

In our analysis we have chosen to fit a single power law 
to the observed field counts. We used a maximum-likelihood 
technique, based on minimising to determine the slope 
(7) and intercept (Cm) of the power law fit. The differential 
counts per square degree per 0.5 magnitude were found to 
be: 



logio dN = (0.331 ± 0.018)1 - (3.24 ± 0.40) 



logio dN 



{0.377 ± 0.028)1/ - (4.73 ± 0.63) 



(2.1) 
(2.2) 



over the range 19.0-22.0 for / and 19.5-22.0 for V (the mean, 
and upper and lower limits on the counts are shown by the 
solid lines in Fig. The errors in the slope and normalisa- 
tion are composed of two components, errors in the fit and 
frame-to-frame errors caused by Poissonian fiuctuations be- 
tween frames. Our slopes are in reasonable agreement with 
those published previously by DPDMD 7/ = 0.34 ± 0.03 
(7 = 19.0-22.5) and 7v = 0.41 ± 0.01 {V = 20.5-23.0). 
Furthermore, the number counts have been found to rise 
line arly in these pass bands to much fainter magnitudes (see 
e.g. small et al. 1995) and so we feel confident in marginally 



extrapolating our field counts to obtain the same depth as 
used in the cluster analysis. 

At this point we wish to quantify the effect gravitational 
lensing by the clusters will have on the surface density of 



background galaxies seen through the cluster. Both clusters 
are known to be strong gravitational lenses and they have 
the potential to increase or decrease th e observed surfac e 
densities of galaxies seen through them (Broadhurst 1996). 
Two competing effect are present: the light paths to distant 
galaxies are deflected by the cluster, resulting in their im- 
ages being displaced radially outward from the cluster centre 
(effectively lowering the background surface density), mean- 
while faint galaxies are brightened by the lensing amplifi- 
ca tion above our ma gnitude limit. Using the methodology 
of Broadhurst (1996 ), we expect the observed number den- 
sity of galaxies, N', to be related to the number density of 
galaxies in the absence of the cluster. No, by: 



^' _ ,2.5.-1 



(2.3) 



where A is the amplification and s is the slope of the field 
counts. Using a singular isothermal sphere approximation 
to the cluster potential and our V and 7 band field count 
slopes of 0.377 and 0.331, we would expect to find a deficit 
of 2.5 per cent and 7 per cent respectively in the background 
population in each passband for our fieldsize. These effects 
are negligible compared with the field-to-field scatter and so 
we will ignore lensing in the ensuing discussion. 



3 ANALYSIS 

In summary, we have measured the differential galaxy counts 
in our cluster fields and we have discussed the photometric 
zero point offsets applied to the AAT catalogues to ensure 
consistency of photometry between the AAT and INT data 
sets. We have determined power law fits to the background 
field counts and hence we are now in a position to subtract 
the number of background galaxies statistically as a function 
of magnitude from the total number of galaxies in our cluster 
fields, to leave only the cluster members. In the following 
sections we will firstly fit single Schechter functions to our 
cluster luminosity functions as is customary, before going on 
to show that improved fits may be obtained using combined 
Gaussian-)- Schechter functions. 



3.1 IVIeasuring the Luminosity Functions 

We show the field-corrected cluster LFs for the clusters 
in both passbands in Fig. ^. These include corrections for 
incompleteness in the cluster frames. One additional sys- 
tematic effect remained which had to be corrected for: in 
a crowded cluster field there is less open area available in 
which to detect the fainter galaxies than in the correspond- 
ing field sample and so some faint background objects will 
be obscured by brighter galaxies (as opposed to just having 
merged isophotes) . This effect will also occur in the field but 
will be less important. We attempted to correct for this by 
computing the available clear field of view for each magni- 
tude interval by subtracting from the total frame area the 
area covered by brighter galaxies. This technique will slightly 
overestimate the compensation factor required, because it 
assumes that all faint objects which happen to occupy the 
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Figure 3. Number counts for the / (upper) and V (lower) images of the two clusters. The counts are per field of view area for each cluster 
(Table |^) and per 0.5 magnitude bin. The solid points have been corrected for incompleteness, the dotted points are the uncorrected 
counts offset by 0.1 mag for clarity. The equivalent field counts from the AAT f/1 data are overlaid. The three solid lines show the best 
fit to the field counts and the ±lo" errors. The dashed line shows the 50 per cent incompleteness limit. Note the excess in the cluster at 
7 ~ 17 and V ~ 18.5 associated with the onset of the bright cluster population. Stars have been removed from all the counts. 



8 Wilson et al. 



Mj Mj 



-24 -22 -20 -18 -24 -22 -20 -18 




16 18 20 22 16 18 20 

I I 




Figure 4. Field corrected cluster LFs, defined as number of galaxies 
both incompleteness and obscuration by brighter objects (see text), 
and have been offset by 0.1 mag for clarity. Overlaid are the best 
Tables | and |. 



per field per 0.5 mag. The solid points and line have been corrected for 
while the dotted points and line are only corrected for incompleteness 
fit Schechter functions for each cluster. The parameters are given in 
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Figure 5. a and M* error contours. The contours are Xmin ^ '^ which corresponds to a 90 per cent confidence ellipse for normally- 
distributed errors and three free parameters. The solid lines are from a Schechter fit to the incompleteness and obscuration corrected 
points. There are three lines corresponding to the best and extreme (itlcr) values of the power law fits to the field counts. The dotted 
lines arc from a Schechter fit to the incompleteness only corrected points. The bold line is from a Schechter fit to the incompleteness and 
obscuration corrected inner/outer test points. 
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same region of sky as a brighter counterpart will be ob- 
scured i.e. that the bright galaxies are "optically-thick" at 
their isophotal diameter. The "true" obscuration correction 
will depend on the surface brightness distributions in the 
bright and faint galaxies, the internal extinction distribu- 
tion in the bright galaxies and the properties of the detection 
algorithm. We claim, therefore, that the two limits we use 
(no correction and the optically-thick case) span the likely 
range. 

As for the field counts, we again used a maximum- 
likelihood method based on chi-squared to measure the a, 
M* and N* parameters in the Schechter function, 



iV(M) dM = fc7v*e'=("+^»"*-")e-=''"*""' 



dM 



where fc = | x (In 10). 

The generalised chi-squared statistic used was 

2 



2 ^Sr^ f logioA^o(i) -logio NE(i 



(3.1) 



(3.2) 



where A''_B(i) is the number of galaxies expected from the 
Schechter function in the i'th bin, and No{i) is the number 
observed. cr(i) is the error for each bin, namely, the Poisson 
error in the cluster-plus-field counts and the Poisson error in 
the field counts, added in quadrature. In addition, we have 
chosen to double the error for the final bin to reflect the 
greater uncertainties associated with final point measure- 
ments and correction values. 

As a simple estimate of the scatter in the Schechter 
fit due to uncertainties in the field counts, we fit for our 
Schechter parameters three times, subtracting firstly the 
field counts derived from our best fit field slope and normali- 
sation parameters , ari d then the extreme (±lcr values) (from 
equations 2.1 and 2.2). This gives three best fit values for the 
three cases. As the Schechter function involves three param- 
eters, the confidence contours are three-dimensional shapes. 
The dotted lines in Fig. ^ show one slice through each 3-D 
contour shape (A''* is held constant at its most likely value). 
Ellipses corresponding to Xmiii+6.3 (the 90 per cent confi- 
dence level for normally distributed errors with three de- 
grees of freedom) error contour are marked. The most likel 
values of the Schechter parameters are shown in Table 
The function itself is shown in Fig. ^ (solid for both incom 
pleteness and obscuration corrected versions and dotted for 
incompleteness corrected only) . To convert from apparent to 
absolute magnitudes we used M — m — Slogj^g D — 25 — K 
where D is the luminosity distance to the cluster {Hq = 50 
kms~^Mpc~^) and K is the K-correction for E/SO galax- 
ies (calculated by convolving our filter responses with the 
spectral energy distribution of a present-day elliptical galaxy 
redshifted to z = 0.18). The combined K-correction and dis- 
tance modulus, m — Af, is 40.50 in I and 40.82 in V. 



3.2 Independent Comparison Check 

We next determined a rough estimate of the cluster LF 
which was entirely independent of the relative magnitude 
scales of the cluster and field data. Here we took advantage 
of the large field of view available in our clusters and the 
strongly peaked distribution expected for galaxies bound to 



the cluster potential. We split our cluster frames into two 
independent radial bins, each covering the same chip area. 
For the various images this translated into radii of ~ 0-0.9 
Mpc and ~ 0.9-1.3 Mpc. We then simply subtracted the dif- 
ferential counts in the outer region from those in the inner 
one. This would completely remove any galaxy population 
which has a flat distribution across the frame, i.e. the fleld 
population, and leave only the peaked cluster galaxies. In 
the absence of luminosity segregation within the clusters, 
the resulting magnitude distribution of galaxies would rep- 
resent the global LF of the cluster]^ 

We show the results of this analysis in Fig. ^ The solid 
points and line have been corrected for both incompleteness 
and obscuration by brighter objects, while the dotted points 
and line have only been corrected for incompleteness and are 
offset by 0.1 mag for clarity. Overlaid on Fig. ^ is the best 
flt Schechter function. The parameters are given in Table ^ 
Note that we might expect to obtain a smaller value for 
N* , the normalisation, using this method, compared to the 
value obtained in the fleld subtraction case. This is because 
some fraction of cluster galaxies will be located in the outer 
annulus and will be subtracted from the total along with 
the fleld galaxies. We see by comparing values of N* from 
Tables ^ and ^ that this is indeed the case. 

The bold lines in Fig. ^ show one slice through each 
error flt to a Schechter function for the incompleteness and 
area corrected points. The errors are much larger for this 
test than the fleld-subtraction one, but we reiterate that as a 
differential test this is not affected by a zero-point magnitude 
error and is sensitive to any population of clustered galaxies 
centred on the cluster centre. 



4 RESULTS AND DISCUSSION 

In Fig. ^ we showed the best flt Schechter functions includ- 
ing both incompleteness and area corrections. In Fig. ^ we 
showed the results of an independent, if less sensitive test. In 
Fig. ^ we made a quantitative comparison of the various flts 
by displaying the error contours associated with each of the 
distributions. Comparing the A665 and A1689 V and / er- 
ror distributions in Fig. M with each other it is clear that the 
best flt parameters for a single Schechter function are only 
in very rough agreement, both between the clusters and for 
the various flts within each cluster. This is surprising given 
the apparent similarity of the LFs in Fig. |^. The reason 
for this is apparent when we compare the best flt Schechter 
functions to the data, especially the V data. Both clusters 
exhibit a faint component which rises very rapidly at faint 
limits (The coarser binning and lower signal to noise in the 
inner/outer comparison smears out this component and re- 
sults in a generally steeper faint end slope to the flts) . A sin- 
gle Schechter function is incapable of simultaneously fltting 
both this feature and the bright end of the LF, as is shown 



t We have tested for luminosity segregation (Lobo et al. 1996) 



by combining the LFs in the inner and outer regions of the two 
clusters and find only weak evidence for a steepening of the faint 
end slope a with radius in the clusters. 
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Figure 6. Estimates of the cluster LPs obtained from differential counts between the central and outer regions of each cluster. The solid 
points and line have been corrected for both incompleteness and obscuration by brighter objects (see text). The dotted points and line 
have been corrected for incompleteness only and offset by 0.1 mag for clarity. Overlaid are the best fit Schechter functions for each case. 
The parameters for the incompleteness and obscuration corrected fits are given in Table 0. 
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Figure 7. Incompleteness and obscuration corrected points for A665 (solid squares) and A1689 (solid triangles). Also shown (open 
squares) is the data from Driver et al (1994) which has been scaled as described in the text. Overlaid to our data is the best fit 
Gaussian+Schechter function (See Table ^ for the parameters). The A1689 data and line has been offset by 0.1 mag. Note the good 
agreement between the form of the LF in the three clusters. 



by the reduced /'^ values for the fits hsted in Table ^ 
These range between 2.1 and 6.6 and indicate that a single 
Schechter function is a poor analytic fit to our data points. 
The disagreement between the fits for the two clusters thus 
arises from the combination of the incorrect functional form 
for the LF and the different contributions of the faint compo- 
nent to the fit due to the slightly different magnitude limits 
of the datasets. 

In order to decide upon a better functional form to de- 



scribe t ]he observed LFs we again plot our V and I data in 
Fig. The A1689 data has been displaced (solid triangles) 
by 0.1 mag from the A665 data (solid squares). However, it 
was not necessary to scale the data vertically (i.e. both clus- 
ters are of approximately equal richness). DPDMD's A963 
R data is also shown displaced and overlaid (open squares). 
We have applied a combined {V — R) colour and distance 
modulus displacement of +0.51 for E/SO galaxies {{R — I) 
of 1.08), and also an empirical normalisation constant. The 
very close similarity between the shapes of the three func- 
tions is immediately apparent, with all three functions ex- 



hibiting a Schechter function form at the bright end and then 
evidence tor an upturn at magnitudes fainter than V ~ 22 
(7 ~ 20). In the following discussion we will focus on the V 
data where this upturn is more prominent. 

In the paper of DPDMD it was also noted that the up- 
turn at the faint end ensured that a single Schechter function 
was not a satisfactory fit to the data. DPDMD decided to fit 
a double Schechter function to their data, fixing the slope 



of the first function (with the brighter knee) to be equal 
to —1.0. Here, instead, we elected to fit a combination of a 
Gaussian function at the bright end of the form: 



N{M) = Ke 



-{M-AI'^""')2/2n 



(4.1) 



and a Schechter function given by equation 3.1 at the faint 
end. We prefer this description over the double Schechter 
function as it has been used e xtensively in i n vestigations of 



the LFs of local clusters (e.g. Metcalfe 1983|; Seeker & Har- 



ris 1996 



Thompson & Gregory 1993; [Biviano et al. 1995) 



including morphological studies (e.g. [Binggeli et al. 198! 



which have shown that giant galaxies tend to follow a Gaus- 
sian distribution and dwarfs a Schechter function. 

As there is some degeneracy between the various pa- 
rameters we have elected to fix the the peak of the Gaus- 
sian, K, to be 10 per cent of the normalisation, A^*, of the 
Schechter function (reasonable values lie in the range 0.08 
- 0.14). We also chose to fix the standard deviation of the 
Gaussian to be a = 1.0, in keep i ng with values determined 
locally for Com a (Metcalfe 1983; Biviano et al. 199^ ; ^eckei 
& Harris 1996). Table H gives the best fit parameters for the 



Gaussian+Schechter function description of the individual 
cluster LFs, as well as those for the samples without the ob- 
scuration correction. As is apparent, this function provides a 
much better description of the data than a single Schechter 
function. This is supported by the reduced values for 
the Gaussian-|-Schechter fits which are all significantly bet- 
ter than the sing le Schecter fits (Table §) . The reduced 
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values are still higher than might be expected, although 
we have not allowed for the eflFects of clustering and have 
low numbers of objects in our bright bins. Nevertheless, the 
Gaussian+Schechter fits allows us to make a robust quanti- 
tative comparison of the form of the LFs in the two clusters. 

We see from Table ^ that there is remarkable consis- 
tency between the M'^'"^*' and M* values for each bandpass. 
In the V band, M""''* is -20.27 for A665 and -20.19 for 
A1689, and in V it is -22.14 for A665 and -22.06 for A1689. 
There is also close agreement between the values of AI* , with 
-18.13 for A665 and -18.50 for A1689 in V, and -20.49 for 
A665 and —20.47 for A1689 in I. In view of this good agree- 
ment between the two clusters we therefore combined their 
LFs and show the fit to this combined dataset in Fig. ^. The 
parameters used may be found in Table |^. 

The main point of interest in comparing the cluster lu- 
minosity functions in V and I is the much steeper faint 
end slope in V compared to I. This difference implies that 
fainter cluster galaxies must have colours which are signifi- 
cantly bluer than the bright ellipticals. In order to quantify 
this effect we plotted the cumulative number of galaxies (us- 
ing our best Gaussian-|-Schechter fits) in each passband as 
a function of magnitude. If both passbands were detecting 
the same objects then comparing the magnitude limits in 
V and / as a function of cumulative density would give the 
typical colours as a function of magnitude in the population. 
We found that the mean (V — /) shifts from {V — I) ^ 1.6 
at / ~ 16 (the colour of the bright cluster ellipticals) to 
(V - /) ~ 0.6 by / = 22 (M = -18.5). The steady blueing 
of the cluster elliptical population as a function of mag- 
nitude, presumably due to metallicity efi^ects, amounted to 
only d{V—I) /dl — 0.07 per mag, and so can only account for 
half of the colour shift. Thus the difference in the faint end 
slopes in the V and I passbands apparently refiects a rapid 
blueing trend in the faint galaxy population in the clusters, 
with the faintest cluster members having the colours and 
luminosities of typical dirr galaxies. 

As we discussed earlier, the Gaussian-|-Schechter pa- 
rameterisation has been used to describe the LF of the Coma 
cluster (an Abell richness class 2 cluster). Here, we com- 
pare the properties of the bright galaxy populations in our 
distant clusters, as described by the Gaussian component, 
with those seen in Coma. The value of M'^'^"* for our com- 
bined cluster sample is M'"'"* = —20.20 in V, very close t o 
that foijind tor Coma: 
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Figure 8. Our simple model for LF evolution in clusters. The 
composite LF of the cluster is shown by the solid line (except 
for the 2 = case where we show the composite dwarf LF), 
the various sub-components are shown by dotted/dashed lines. 
Overlaid on these are the observed data points from Fig. ^ We 
show in the bottom panel the dwarf elliptical LF derived for Virgo 
by Impey et al. (1988) converted into V using (B — V) = 0.6 
(applicable to dEs) , where we have increased the size of the error 
bars at the faint end to take into account the uncertainties in the 
corrections applied by Impey et al. 



-20.20"t piviano et al. 199bD and — 2U.lb ( ]Seckor k Harris 
1996 |) ([ converting all measurements to the V band using 
{B — V^ = 1.0). The scatter amongst the various analyses of 
Coma indicates that we ought not to put too much weight 
on this comparison, but it is interesting to note that the 
distant clusters appear to be ~ 0.1 brighter in rest-frame V 
compared to the local value. This degree of brightening is 
not unreasonable out to z = 0.18 in a population of p assively 
evolvin g el liptical galaxies formed at high redshift (Barger 
et al. iQ. 

Turning to the faint cluster population, parameterised 
by the Schechter function, we find that locally the values 
of the characteristic luminosit y Al* are more poorly defined 
and have larger errors; —18.1 ( Biviano et al ), —18.7 ( Becker 
& Harrj^) and —18.0 for the dE population in the Virgo clus- 



20.03 ([Thompson & Gregory 1993|), Harris| ) and possibly even steeper (a 



ter (Binggeli et al. 1985). Again these encompass the value 
observed in our distant clusters, Al* — —18.5, indicating a 
similar lack of strong evolution in the characteristic luminos- 
ity of this population to that seen for the brighter galaxies. 
However, when we compare the faint end slope, q, for the 
faint cluster population we find a large difference between 
the distant and local clusters. The faint end s lope in Coma 

lies in the range ~ —1.1 1.4 (Biviano et al. and Becker & 

1.7 for the dE pop- 
ulation of Virgo (Impey et al. 1985) ), compared to a ~ —2 
for our distant clusters Q 

The consistency of the upturn at the faint end in the 
three clusters available to us, particularly in the V band, 
suggests that this feature is common to many distant rich 
clusters. The fact that the upturn appears much less dra- 
matic in the I band indicates that faint end population is 
significantly bluer than the typical bright cluster galaxies. 
To obtain an overview of the implications of these observa- 
tions we have constructed a simple model which takes the 

t The divergent slope found in our fits can easily be reduced 
by selecting a different characteristic luminosity in the Schechter 
part of the LF. 
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broad features of the observed z ^ 0.2 V band LF discussed 
above and predicts both the z ~ 0.2 / band LF and the 
evolved V band LF which we would expect to observe lo- 
cally. The key feature of this model is the division of the 
faint component of the cluster LF into two equal popula- 
tions, dwarf ellipticals (dE) and dwarf irregulars (dirr). The 
LFs for both of these groups follow a Schechter form, but 
the different star formation histories expected for the two 
populations result in significantly different evolution to the 
present day. 

In the top panel of Fig. ^ we replot our V band data 
from Fig. ^ We also show the constituent LFs of the three 
populations assumed in our simple model, and the resultant 
LF. The giant elliptical galaxies (gE) have been assigned the 
same M'^'^"* and A^* values as discussed in Table |^. The dwarf 
elliptical and dwarf irregular populations have equal normal- 
isations and characteristic luminosities {My ~ —17.8) but 
have radically different faint end slopes. The dwarf ellipti- 
cals have been assigned a value of a ~ —1.2 (chosen because 
this was the value observed for our I band LF). The dwarf 
irregular population has a much more steeply rising slope 
(with a ~ —2). 

We next used the colours expected for these three pop- 
ulations to predict the expected form of the z ^ 0.2 / 
band LF. Wo uaod tho oboorvod coloura of the giant and 



dwarf allipticals taken from fits to the E/SO sequence in 
the clusters and estimated at their characteristic magnitudes 
{{V - J) ~ 1.5 for Es, and (F - /) ~ 1.2 for dEs) to trans- 
form their LFs from V to /. For the dIrr population we 
adopted the colour of a dIrr galaxy (close to fiat spectrum), 
{V — I) ~ 0.6. The three transformed LFs are shown in 
Fig. ^ as is the expected composite I band LF. As can be 
seen this is in reasonable agreement with the observations 
which we have overplotted. 

Finally we used our model to predict what present day 

V band LF we might expect to observe in a rich cluster. 
The two faint end populations will evolve very differently 
with time. The dwarf ellipticals will evolve passively in a 
similar manner to the bright cluster ellipticals (with fading 
by about 0.1 mags in rest-frame V between z — 0.2 and 
today). The very blue (star forming) dwarf irregulars, how- 
ever, are assumed to have their star formation stopped soon 
after we observe them (around ~ 3 Gyrs ago), either due 
to a violent truncation mechanism (such as ram-pressure 
stripping of their disks or a galactic wind) or by the slower 
removal of the gas-reservoirs in their halos. This population 
will then fade rapidly (by upwards of 2-3 mags in rest -frame 

V ( Barger et al. 1996 ), and possibly up to ~ 5mag (Babul 
& Ferg fson 199'^ An even more drastic possibility would 
be that the galaxies are partially or totally disrupted due to 
tidal interactions in the cluster and either become low sur- 
face brightness dwarfs (the bl ue and hence relativ ely young 
objects discovered in Virgo by Impey et al. (1988) or be re- 



moved from the scene altogether (c.f. Babul & Rees 1992) 



Whatever the mechanism, the dwarf irregulars are assumed 
to have faded dramatically (3 mags) by the present day. We 
show the individual LFs in the bottom panel of Fig. M and 
also the composite dwarf LF (the fade d dIrr-l-dE po pula- 
tion). We have also plotted in Fig. ^ the [mpey et al. LF of 



equately described using our combined dwarf population. It 
is interesting to note also that the region around My ~ — 16 
where the dIrr beg in to appear i n reasonable numbers is also 
the point at which Impey et al. report the breakdown in the 



surface-brightness-luminosity correlation in their data, indi- 
cating that these fainter objects probably arise from a very 
mixed group of progenitors. 

In summary, the significant difference in the colours of 
the dwarf ellipticals and dwarf irregulars implies that the 
former dominate in the / band LF of the distant clusters to 
the limit of our data (although not much fainter) and the 
latter determine the faint end in the V band. The differential 
evolution in the dwarf irregulars results in them fading far 
down the LF by the present day - leaving the intermediate 
regime My ~ — 18 dominated by the dwarf ellipticals. Below 
Mv ~ — 16 remnants of the dwarf irregulars begin to be seen 
in large numbers resulting in an upturn in the LF. 

Finally, we used our model to predict the early-type 
dwarf to giant gal axy ratio (EDGR) in o ur distant clusters. 
Local studies i.e. Seeker & Harris (1996) find an EDGR of 
1.8 ± 0.6 in the Coma cluster to a limit of Mb ~ —15.5 
(My ~ —16.1), similar to t he value of 2.1 found for Virgo by 



Ferguson & Sandage (1991). Note that we co mpared our re- 



sul ts o nly with the brightest subsample from Becker & Har- 
ris] and Ferguson fc Sandage because this required the least 



extrapolation beyond the limit of our data. If we calculate 
the EDGR from the dE and gE populations in our model at 
z ~ 0.2, we obtain an EDGR of 2.3 for {K/N*) =0.1 (rang- 
ing from 3.3 for {K/N*) = 0.08 to 1.9 for [K/N*] = 0.14). 
Little evolution is expected in our model for the EDGR 
brightward of My ~ —16 as the faded remnants of the 
dwarf irregulars do not appear brighter than this at the 
present day (It is not clear whether the proposed popula- 
tion of dwarf irregulars could metamorphose into e.g. dwarf 
spheroidals in the interim between z ~ 0.2 and the present 
day. The point we stress here is that with our assumption of 
3 mags fading the remnants cannot affect the EDGR bright- 
ward of My ~ — 16 whatever form they take). Thus from 
our model we expect an EDGR of 2.3lJ'4 local clusters, 
consistent with the values observed. We conclude that there 
is currently no convincing evidence for strong evolution in 
EDGR out to z - 0.2. 

Finally, turnin g to the "galaxy harassment" model of 
Moore et al. (199^ we see that the predicted evolution of 
the EDGR will depend critically upon the relative effect 
of disk stripping and central sta rbursts on th e luminosity 
of the remnant dwarf spheroidal. 



Moore et al 



dwarf ellipticals in the Virgo cluster. Glearly, this can be ad- 



suggest that 

the infalling field galaxies are somewhat brighter than the 
dIrr population in our models {My ~ —19). If the stellar 
luminosity lost due to stripping is less than a few times 
that produced by the starburst, the remnants would appear 
brighter than My ~ —16 in the local clusters and hence 
the present day EDGR would be higher than that seen at 
intermediate redshift, in contradiction to the observations. 
Glearly there is considerable scope for obtaining agreement 
between the harassment model and the observations, but we 
are encouraged by the future possibilities of closely compar- 
ing theoretical models with observations of galaxy evolution 
in rich clusters. 
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5 CONCLUSIONS 

(i) We have presented deep V and / photometry of A665 
and A1689, two rich, X-ray luminous clusters at z — 0.18. 
This photometry reached a limiting apparent magnitude of 
/ = 22.5, equivalent to an absolute magnitude of 7 = —18.0, 
or M* + 5, at the cluster redshifts. We analysed the data to 
provide differential number counts of ~ 1500-2000 galaxies 
in the fields of each of these two clusters. 

(ii) We used an independent V and I field survey reaching 
similar depths to remove the field contamination from our 
cluster fields. We undertook extensive observations to cali- 
brate the photometry between the field and cluster datasets. 
This gave us confidence in our field subtraction, and in the 
resulting cluster LF. 

(iii) The LFs we derived in our two clusters also showed 
remarkable consistency with each other. Both showed knees, 
a flat region and then an upturn at the faint end (beginning 
around My = —19). The upturn in the cluster LFs was 
more marked in V than / (a ~ —2 versus a ~ —1.2) indi- 
cating that the faint population rapidly becomes bluer with 
decreasing luminosity. By means of a simple test using cumu- 
lative plots of the galaxies in each band we showed that the 
mean colour was (V — I) ~ 0.6 by / ~ 22 in contrast with the 
mean colour of the bright spheroid population, (V — I) ~ 1.6. 
Moreover, the form of the LF in our two clusters was in very 
go od agreement with that found for a single distant cluster 
by [Driver et al. (1994 . 

(iv) We constructed a simple model which described the 
gross features of the LFs observed in our distant clusters. 
The main feature of this model was the partitioning of the 
faint cluster component into a quiescent red dwarf elliptical 
population and a star forming blue dwarf irregular popu- 
lation. Evolving these populations forward to the present 
day we could fit the form of the dwarf luminosity function 
observed in Virgo if we assumed substantial, but not unrea- 
sonable, fading in the dirr population. 

(v) If the cessation of star formation and subsequent fad- 
ing/disruption of these blue dwarf galaxies is associated with 
extensive gas loss from these systems, then the population 
seen in our distant clusters has a sufficiently steep faint end 
slope (q ~ — 2) to be the source of all the X-ray gas seen in 
local rich clusters (Trentham 1994). 



(vi) Extending observational studies such as the one pre- 
sented here to higher redshifts will provide further con- 
straints on the changing form of the giant and dwarf galaxy 
populations in clusters, and hence the consequences for the 
evolution of the X-ray emission of the clusters. 
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Table 1. INT cluster observations 



Target 


Filter 


Chip 


Scale 


Field 


Texp 


FWHM 


miim 


m(1ct) 


Reddening 


Ngal 








C'/pix) 




(ksec) 


(") 


(50%) 


(mag/D") 




(< mso) 


A665 


I 


FORD 


0.37 


12.34' X 12.30' 


20.5 


1.7 


22.0 


25.3 


0.007 


1689 


A1689 


I 


EEV 


0.54 


11.29' X 10.61' 


9.8 


1.8 


22.5 


26.1 


0.063 


2035 


A665 


V 


FORD 


0.37 


12.34' X 12.30' 


17.0 


2.0 


23.5 


26.7 


0.015 


1659 


A1689 


V 


FORD 


0.37 


11.29' X 10.61' 


18.0 


2.1 


23.5 


26.9 


0.098 


1615 



Table 2. SExtractor detection parameters 



Cluster 


Filter 


Conv. 


Min. Area 
(pixels) 


Mthrcsh 

(mag/n") 


A665 


I 


3x3 


6 


24.5 


A1689 


I 


3x3 


6 


25.1 


A665 


V 


3x3 


6 


26.3 


A1689 


V 


1x1 


6 


26.3 



Turner J. A., Phillipps S., Davics J. L, Disney M. J., 1993, MN- 

RAS, 261, 39 
Tyson J. A., Fischer P., 1995, ApJ, 446, L55 
Tyson J. A., Valdes F., Wenk R., 1990, ApJ, 349, LI 
Valdes F., Tyson J. A., Jarvis J. F., 1983, ApJ, 271, 431 
Whitmore B. C, Gilmore D. M., Jones C, 1993, ApJ, 407, 489 

This paper has been produced using the Royal Astronomical 
Society/Blackwell Science MJ^X style file. 
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Table 3. Schechter function parameter fits for incompleteness corrected field subtracted points. The columns show the cluster, filter, 
faint end slope, apparent magnitude at the knee, absolute magnitude at the knee, normalisation, , and reduced values. The best 
fit parameters found by subtracting the mean value (first row) and then the ±la error values of background counts (second and third 
rows) are also shown. 



Cluster 


Filter 


a 


m* 


M* 


N* 






A665 


/ 


-1.065 


18.11 


-22.39 


127.15 


24.67 


2.74 






-1.083 


18.04 


-22.46 


143.03 


24.54 








-1.060 


18.18 


-22.32 


103.57 


25.41 




A1689 


/ 


-1.321 


17.42 


-23.08 


62.63 


21.34 


2.13 






-1.303 


17.40 


-23.10 


73.40 


22.40 








-1.366 


17.40 


-23.10 


46.76 


20.54 




A665 


V 


-0.843 


19.64 


-21.18 


240.54 


49.36 


5.49 






-0.709 


19.73 


-21.09 


296.56 


40.87 








-1.079 


19.46 


-21.36 


154.46 


57.58 




A1689 


V 


-1.500 


18.83 


-21.99 


57.25 


33.33 


3.70 






-1.382 


19.00 


-21.82 


85.73 


30.48 








-1.647 


18.56 


-22.26 


30.93 


34.92 





Table 4. Schechter function parameter fits for incompleteness and obscuration corrected field subtracted points. The column headings 
are as in Table H. 



Cluster 


Filter 


a 


m* 


M* 


N* 


x' 


xV^ 


A665 


/ 


-1.188 


18.04 


-22.46 


123.00 


26.03 


2.89 






-1.185 


17.98 


-22.52 


139.02 


25.87 








-1.222 


18.04 


-22.46 


95.77 


26.69 




A1689 


/ 


-1.430 


17.30 


-23.20 


57.04 


22.12 


2.12 






-1.406 


17.28 


-23.22 


66.52 


23.01 








-1.466 


17.37 


-23.13 


47.38 


21.42 




A665 


V 


-1.167 


19.37 


-21.45 


178.36 


59.17 


6.57 






-1.042 


19.43 


-21.39 


233.79 


52.01 








-1.370 


19.00 


-21.82 


96.92 


64.57 




A1689 


V 


-1.607 


18.68 


-22.14 


51.40 


35.96 


4.00 






-1.514 


18.84 


-21.98 


74.56 


34.07 








-1.715 


18.48 


-22.34 


31.70 


36.80 





Table 5. Schechter function parameter fits for incompleteness and obscuration corrected radial inner/outer test points. The column 
headings are as in Table H. 



Cluster 


Filter 


a 


m* 


M* 


N* 


x' 




A665 


/ 


-1.235 


17.49 


-23.01 


40.96 


1.27 


0.14 


A1689 


/ 


-1.034 


17.43 


-23.07 


139.84 


10.24 


1.02 


A665 


V 


-1.264 


18.82 


-22.00 


58.07 


4.42 


0.49 


A1689 


V 


-1.681 


18.09 


-22.73 


19.10 


5.67 


0.63 
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Table 6. Gaussian+Schechter function parameter fits. The top table gives the fits for the incompleteness and obscuration corrected 
catalogues, while the lower is for the incompleteness corrected case only. The columns show the cluster, filter, apparent magnitude at the 
peak of Gaussian function, absolute magnitude at the peaJc of Gaussian function, normalisation of Schechter function, faint end slope 
of Schechter function, apparent magnitude at the knee of Schechter function, absolute magnitude at the knee of Schechter function, i 
and reduced values. Note that a, the standard deviation of the Gaussian function is held constant at 1.0 and K, the normalisation 
of the Gaussian function is required to equal 0.1 X N* . 



Cluster 


Filter 


^ccnt 


^ccnt 


N* 


a 


m* 


M* 






A665 


/ 


18.36 


-22.14 


312.21 


-0.991 


20.01 


-20.49 


12.01 


2.00 


A1689 


/ 


18.44 


-22.06 


319.99 


-1.180 


20.03 


-20.47 


11.89 


1.49 


Combined 


I 


18.66 


-21.84 


343.98 


-1.160 


20.23 


-20.27 


12.32 


1.76 


A665 


V 


20.55 


-20.27 


823.73 


-2.295 


22.69 


-18.13 


9.51 


1.36 


A1689 


V 


20.63 


-20.19 


689.95 


-1.870 


22.32 


-18.50 


9.99 


1.43 


Combined 


V 


20.62 


-20.20 


776.95 


-2.092 


22.55 


-18.27 


14.68 


2.10 


Cluster 


Filter 


^ccnt 


^ccnt 


JV* 


c\ 


m* 


M* 


X^ 


xV'^ 


A665 


/ 


18.22 


-22.28 


268.05 


-0.841 


19.97 


-20.53 


11.68 


1.95 


A1689 


/ 


18.42 


-22.08 


290.07 


-1.025 


20.18 


-20.32 


12.01 


1.50 


Combined 


/ 


18.64 


-21.86 


309.18 


-1.017 


20.36 


-20.23 


12.31 


1.76 


A665 


V 


20.48 


-20.34 


673.15 


-2.464 


22.81 


-18.01 


9.17 


1.31 


A1689 


V 


20.60 


-20.22 


572.42 


-1.956 


22.45 


-18.37 


8.92 


1.27 


Combined 


V 


20.56 


-20.26 


635.45 


-2.230 


22.66 


-18.16 


12.69 


1.81 



